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Summary 

Although  stringers  are  used  primarily  as  stiffen¬ 
ers,  they  also  can  make  damaged  structures  fail-safe 
or  damage  tolerant.  Assessment  of  the  damage  toler¬ 
ance  of  structures  weakened  by  cracks  is  aided  by 
knowledge  of  stress-intensity  factors.  In  this  pa¬ 
per,  the  stress-intensity  factor  is  determined  for  a 
cracked  orthotropic  sheet  adhesively  bonded  to  an 
orthotropic  stringer.  Since  the  stringer  is  modeled  as 
a  semi- infinite  sheet,  the  solution  is  most  appropriate 
for  a  crack  tip  located  near  a  stringer  edge.  Both  ad- 
herends  are  treated  as  homogeneous,  orthotropic  me¬ 
dia  which  are  representative  of  many  fiber-reinforced 
composite  materials.  It  is  assumed  that  the  ad- 
herends  are  in  a  state  of  plane  stress  and  the  adhesive 
is  in  pure  shear.  By  use  of  Green’s  functions  and  the 
complex  variable  theory  of  orthotropic  elasticity  de¬ 
veloped  by  Lekhnitskii,  a  set  of  integral  equations  is 
obtained.  The  integral  equations  are  replaced  by  an 
equivalent  set  of  algebraic  equations,  which  is  solved 
to  obtain  the  shear-stress  distribution  in  the  adhesive 
layer.  From  these  adhesive  stresses,  equations  for  the 
stress- intensity  factors  at  both  crack  tips  are  found. 

A  parametric  study  is  also  conducted  to  deter¬ 
mine  the  sensitivity  of  the  system  to  material  prop¬ 
erties  and  specimen  configuration.  The  parameter 
having  the  greatest  influence  on  the  stress- intensity 
factors  is  the  distance  from  the  near  crack  tip  to  the 
edge  of  the  stringer.  Unless  the  crack  tip  is  very  close 
to  or  under  the  stringer,  the  stress- intensity  factor  is 
approximately  that  of  an  unstiffened  sheet.  How¬ 
ever,  as  the  crack  propagates  beneath  the  stringer, 
the  stress-intensity  factor  decreases  significantly.  In¬ 
creasing  the  stringer  stiffness  or  the  adhesive  stiff¬ 
ness  also  results  in  a  decrease  in  the  stress-intensity 
factor. 

Introduction 

Because  of  their  high  strength  and  light  weight, 
composite  materials  are  finding  increasing  applica¬ 
tion  in  aerospace  structures.  In  general,  the  struc¬ 
tural  configurations  used  in  composite  airplanes  have 
been  very  similar  to  the  sheet-stringer  construction 
currently  used  in  metal  airplanes.  In  metal  airplanes, 
stringers  are  effective  in  making  damaged  structures 
fail-safe  or  damage  tolerant.  For  this  reason,  the 
interaction  of  a  cracked  sheet  and  a  stringer  is  an 
important  problem  which  has  been  investigated  by 
many  authors. 

Romualdi,  Frasier,  and  Irwin  (ref.  1)  considered 
two  problems:  a  crack  centrally  located  under  a  riv¬ 
eted  stringer  and  a  crack  centrally  located  between 
two  riveted  stringers.  In  their  formulation  of  the 
problem,  only  the  two  rivets  nearest  the  crack  were 


considered  effective,  Sanders  (ref.  2)  considered  the 
problem  of  a  crack  positioned  symmetrically  beneath 
a  continuously  attached  stringer.  He  employed  a 
shear-lag  line  stringer  and  assumed  the  sheet  to  be  in- 
extensible  in  the  direction  parallel  to  the  crack  to  ob¬ 
tain  the  stringer  stress-concentration  factor  and  the 
crack-tip  stress-intensity  factor.  Greif  and  Sanders 
(ref.  3)  removed  Sanders’  previous  inextensionality 
assumption  and  symmetry  requirements.  They  con¬ 
cluded  that  the  stringer- induced  reduction  in  the 
stress-intensity  factor  was  a  localized  effect  and  that 
the  stringer  stress-concentration  factor  was  largely 
insensitive  to  sheet  stiffness  in  the  direction  parallel 
to  the  crack.  Bloom  and  Sanders  (ref.  4)  reported 
an  analysis  of  a  discretely  attached  stringer  which 
may  be  broken.  Poe  (ref.  5)  examined  the  problem 
of  a  cracked  sheet  stiffened  by  evenly  spaced,  riv¬ 
eted  longitudinal  stringers.  The  crack  was  assumed 
to  originate  at  a  rivet  location  and  to  grow  beneath 
intact  stringers.  An  extension  of  that  work  (ref.  6) 
considered  the  influence  of  broken  stringers. 

Arin  (ref.  7)  examined  the  effect  of  the  debond¬ 
ing  of  an  infinitely  long  line  stringer  on  the  stress- 
intensity  factor  at  the  crack  tip.  He  assumed  that 
the  stringer  was  adhesively  bonded  to  an  isotropic 
sheet  along  a  line  perpendicular  to  the  crack.  He 
concluded  that  the  stringer  exerts  little  influence  on 
the  stress-intensity  factors  unless  it  is  quite  close  to 
the  crack  tip  and  the  debond  length  is  less  than  twice 
the  crack  length. 

Norris  (ref.  8),  using  a  complex- variable  approach, 
represented  the  finite-width  stringer  with  an  array  of 
line  stiffeners.  The  stringer  was  divided  into  strips, 
and  each  strip  was  represented  by  a  line  stiffener 
attached  to  the  sheet  at  discrete  points,  with  no 
coupling  between  adjacent  line  elements.  Because  the 
analysis  did  not  model  debonding  of  the  adhesive, 
it  predicted  unrealistically  high  loads  in  the  line 
stiffener  nearest  the  crack  tip. 

Experimental  work  has  shown  that  as  the  crack 
tip  approaches  the  stringer,  debonding  can  start  and 
propagate  through  the  adhesive.  Although  Arin 
(ref.  7)  included  adhesive  debonding  in  his  analysis, 
he  modeled  the  stringer  as  a  line,  and  thus  was 
unable  to  account  for  partial  debonding  across  the 
width  of  the  stringer.  Anderson,  Hsu,  and  McGee 
(ref.  9)  modeled  partial  debonding  in  their  analysis 
of  a  cracked  sheet  reinforced  by  a  bonded  doubler. 
However,  they  were  using  a  finite- element  analysis, 
which  would  be  very  cumbersome  and  expensive  to 
use  when  varying  parameters  such  as  crack  length 
and  debond  area. 

A  realistic  analysis  is  needed  that  examines  the 
effect  of  growing  debond  area  on  crack  propagation 
in  adhesively  bonded  structures.  The  purpose  of  this 


paper  is  to  formulate  such  an  analysis  and  to  use  it 
to  study  orthotropic  sheets  reinforced  by  orthotropic 
stringers.  To  meet  this  objective,  the  problem  has 
been  divided  into  two  parts.  In  the  first  part,  the 
problem  is  formulated  with  the  assumption  of  a  linear 
elastic  adhesive  with  no  debonding,  and  the  stringer 
is  modeled  as  a  semi-infinite  sheet  bonded  to  an 
infinite  sheet  containing  a  crack.  In  the  second  part, 
the  effects  of  adhesive  nonlinearity  and  debonding  are 
included.  Although  both  parts  have  been  solved,  for 
brevity  only  the  first  part  of  the  problem  is  presented 
here. 

Nomenclature 


a 

half-length  of  crack,  m 

-\-Ci 

right  crack  tip,  m 

—a 

left  crack  tip,  m 

Aj,Bj,  Cj,Dj 

constant  coefficients  in  layer  j 
stress  functions,  N/m^ 

b 

distance  from  edge  of  stringer 
to  center  of  crack,  m 

^jk 

complex  constants  (j,  fc  —  1,  2) 

D 

domain  of  integration 

Ex,  Ey 

Young’s  moduli  in  the  x-  and 
y-directions,  respectively.  Pa 

FD{w,z,a) 

complex-variable  function 
used  in  displacement-potential 
functions 

FF{w,a) 

complex-variable  function  used 
in  calculation  of  left  crack-tip 
stress- intensity  factor 

FG{w,z,a) 

complex-variable  function  used 
in  displacement  functions  of 
layer  2 

FH{z,w,a) 

complex- variable  function  used 
in  stress  potentials 

Fi[w,a) 

complex- variable  function  used 
in  calculation  of  right  crack-tip 
stress-intensity  factor 

/oi,/o2 

stress  functions,  m^/N 

Gxy 

orthotropic  shear  modulus,  Pa 

Gs 

shear  modulus  of  adhesive 
layer,  Pa 

h, 

thickness  of  layer  j,  m 

i 

=  v^ 

h 

normal  (mode  I)  component  of 
stress-intensity  factor,  Pa\/m 

N 

number  of  collocation  points; 
number  of  cells  in  the  domain 

complex  constants  for  or¬ 
thotropic  sheet  [j  =  1}2), 

1/Pa 

Qjkj^jk 

complex  kernels  used  in 
displacement  functions  for 
layers  1  and  2,  respectively 
(y,fc=l,2),l/Pa 

Sjk 

complex  kernels  used  in 
integral  equations  (j,  fc  =  1,2), 
l/(Pa-m) 

Tj 

complex  kernels  used  in  stress- 
intensity  factors  (/  =  1,  2,  3, 

4) 

X-  and  ^-displacements  of  layer 

displacements  in  layer  1  due  to 
crack-face  pressure,  m 

displacements  in  layer  1  due  to 
concentrated  forces,  m 

(^)  (^) 

displacements  in  layer  1  due  to 
distributed  body  forces,  m 

displacements  in  layer  2  due  to 
a  pair  of  concentrated  forces, 
m 

^2,  V2 

displacements  in  layer  2  due  to 
distributed  body  forces,  with 
respect  to  rotated  coordinate 
system,  m 

Wj 

complex  variable  (/  =  1,  2,  3, 
4),  m 

body  forces  acting  on  layer  j 
[j  =  1,  2),  N/m^ 

layer  2  body  forces  in  rotated 
coordinate  system,  N/m^ 

Cartesian  coordinates,  m 

^o^yo 

coordinates  of  load  point,  m 

rotated  coordinate  system,  m 

^3 

complex  variable  (j  —  1,  2),  m 

^0 

location  of  load  point,  m 

AXj2,  ^yn 

incremental  distances 
(n=  l,2,...,iV),  m 

2 


My 

roots  of  the  characteristic 
equation  (j  —  1,  2) 

^xy  5  ^yx 

Poisson’s  ratios 

pressure  applied  to  crack  face, 
Pa 

'^X')  '^y 

shear  stresses  in  adhesive 
layer,  Pa 

stress  potentials  for  or¬ 
thotropic  materials.  Pa 

4>i{z)AU^) 

displacement  potentials  due 
to  a  pair  of  load  points  in  the 
cracked  sheet,  m 

f2l  (2) 

01  A^))02  displacement  potentials  due 

to  a  pair  of  load  points  in  a 
half-plane,  m 

Superscripts: 

(j),  (fc)  layer  number  {j,k  =  1,2) 

A  bar  over  a  symbol  indicates  a  complex  conju¬ 
gate;  thus,  if  z  =  X  +  iy,  then  z  =  x  —  iy.  A  prime 
on  a  symbol  indicates  reference  to  the  rotated  coor¬ 
dinate  system,  except  in  the  stress  and  displacement 
potentials,  for  which  a  prime  indicates  differentia¬ 
tion. 


Formulation  of  the  Problem 

Consider  the  sheet-stringer  configuration  shown  in  figure  1.  (Hereinafter,  the  semi- infinite  sheet  will  be 
referred  to  as  a  stringer.)  The  stringer  and  the  sheet  are  bonded  together  by  an  adhesive  layer  of  constant 
thickness  /13.  Here  the  bond  is  assumed  to  be  intact.  The  crack  faces  are  subjected  to  a  uniform  pressure  cro 
and  the  stress  state  at  infinity  is  zero,  as  shown  in  part  A  of  figure  2.  The  problem  is  solved  with  the  following 
assumptions: 

1.  The  sheet  (layer  1),  the  stringer  (layer  2),  and  the  adhesive  (layer  3)  are  homogeneous  and  linearly 
elastic. 

2.  The  sheet  and  the  stringer  are  dissimilar,  orthotropic  materials  with  principal  directions  of  orthotropy 
being  oriented  parallel  and  perpendicular  to  the  crack  in  layer  1, 

3.  The  thicknesses  of  the  sheet  and  the  stringer  are  small  compared  with  the  in-plane  dimensions,  so  that 
both  layers  are  considered  to  be  in  a  state  of  plane  stress. 

4.  The  surface  shear  transmitted  through  the  adhesive  is  assumed  to  act  as  a  body  force  on  the  adherends. 

5.  The  thickness  of  the  adhesive  is  small  compared  with  the  thicknesses  of  the  sheet  and  the  stringer;  thus, 
the  adhesive  layer  is  treated  as  a  shear  spring  rather  than  as  an  elastic  continuum. 

Assumption  5  leads  to  the  following  continuity  of  displacement  equations: 


Ul-U2  =  —Tx 


V1-V2  — 


(1) 


where  ui^vi  and  U2,  V2  are  the  x-  and  y-components  of  the  in-plane  displacement  vectors  in  layers  1  and  2, 
and  Ty  are  the  components  of  the  adhesive  shear  stress,  and  /13  and  G3  are  the  adhesive  thickness  and  shear 
modulus,  respectively.  From  assumption  3,  the  two  sets  of  body  forces  (forces  per  unit  volume)  that  act  on 
layers  1  and  2  can  be  written  as  follows: 


Ai-- 


hi 


(2) 


Figure  2  shows  how  the  problem  is  broken  into  its  component  parts.  The  displacements  in  the  sheet,  shown 
in  part  B  of  figure  2,  and  in  the  stringer,  shown  in  part  C  of  figure  2,  are  determined  individually  as  follows. 


Displacements  in  the  Sheet  (Layer  1) 

The  forces  acting  on  the  sheet  (layer  1)  are  shown  in  figure  3.  The  cracked  sheet  is  subjected  to  the  body 
forces  Xi  and  Yi  distributed  over  the  domain  D  and  to  the  crack-face  pressure  ctq.  The  problem  shown  in 
part  A  of  figure  3  is  equivalent  to  the  superposition  of  the  problems  shown  in  parts  B  and  C  of  figure  3. 
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If  and  denote  the  displacements  in  the  sheet  due  to  the  loadings  shown  in  parts  B  and 


C,  respectively,  of  figure  3,  the  total  displacements  in  the  sheet  can  be  written  as 


(?>)  .  (c) 

=  '^1  + 


W  ,  (c) 


(3) 


If  we  follow  Lekhnitskii  (ref.  10)  and  ignore  an  arbitrary  rigid-body  displacement  vector,  the  displacements 
and  due  to  the  crack-face  pressure  (Jq,  shown  in  part  B  of  figure  3,  can  be  found.  The  complete 
expressions  for  and  are  given  in  appendix  A. 

(c)  (c)  .  . 

Determination  of  the  displacements  ^  and  ^  due  to  the  distributed  body  forces  shown  in  part  C  of 
figure  3  starts  with  the  known  displacements  due  to  a  concentrated  force  located  at  the  point  =  ^0  +  ^1/0* 

Since  the  body  forces  Xi  and  Yi  are  continuous  functions  of  zq  in  the  domain  D,  the  displacements  due  to  the 

(c)  fc) 

concentrated  forces  can  be  used  as  Green’s  functions  to  determine  ^  and  .  Thus,  these  displacements 
can  be  written  as  follows: 

u^l\x,y)  =  j j  \Qn  {x,y,XQ,yQ)Xi  (xq,  2/o)  +  <3l2  (a:,  y,  a;o,  2/o)  U  (a;o,yo)]  dx^dyo 

D 

,^-jj  [<321  {x,y,xo,yo)Xi  (a:o,yo)  +  <322  {x,y,XQ,yQ,)Yi  ixo,yo)]  dxodyo 


(c) 


D 


(4) 


where  Qii^Qi2^Q2i^  and  Q22  are  derived  in  appendix  A. 

Displacements  in  the  Stringer  (Layer  2) 

The  only  forces  acting  on  the  semi-infinite  stringer  (layer  2)  are  the  distributed  body  forces  X2and  Y2  shown 
in  part  C  of  figure  2.  In  a  procedure  similar  to  that  used  in  reference  11,  the  displacements  due  to  a  single 
pair  of  body  forces  X2  and  located  at  point  Zq  ~  Xq  Y  iy^  (see  fig.  4),  are  used  to  form  the  displacements 
^2  and  V2  dne  to  the  distributed  loading.  Since  the  body  forces  X2  and  Y2  are  continuous  functions  of  in 
the  domain,  the  displacements  due  to  the  concentrated  force  can  be  used  as  Green’s  functions  for  the  total 
displacements  and  Thus,  these  displacements  can  be  written  as  follows: 

y'^  J  j  [-^11  y'^  ^0’  y'o)  (^O)  y'o)  +  din  {x',  y',  Xq,  yg)  Y2  (xg,  yo)]  dxg  dy'g 

D 

'V2{x',y')  =  j  j  [i?2i  {x'  ,y\x'g,yg)  X'2  (xg,yg)  +  i?22  (x',y',  XQ,yg)  Y2  (xg,yo)]  dxgdy'g 


D 


(5) 


The  expressions  for  i^l25 -^21?  3,nd  R22  are  derived  in  appendix  B. 

The  displacements  U2  and  U2,  required  for  the  loading  shown  in  part  C  of  figure  2,  are  calculated  by  shifting 
and  rotating  the  displacements  ^^2  and  V2  according  to  the  following  relations: 


where 


U2{x,y)  =  -v'2{x’,y')  V2ix,y)  =  1x2(0;',?/') 
y'  =  b  —  X 


(6) 


X  =  y 


X'2  =  >2  =  -X2 

Complete  details  on  the  derivation  of  the  displacements  for  the  stringer  are  presented  in  appendix  B. 

Integral  Equations  for  and  Ty 

Substituting  equations  (3)  and  (6)  into  equations  (1)  yields  the  following  system  of  integral  equations  for 
the  unknown  functions  Tx  and  r^: 


4 


^Tx{x,y)  +  j  j  I  j^^^(a:,y,a:o,yo)  + -^  (a:'iy',a:o,2/o) 

D 


:  (a:o,yo) 


+ 


Qi2 

hi 


{x,  y,  xo,  yo)  -^{x',y',x'Q,  y'o)  Ty  (xq,  j/o)  |  dxQ  dyo  =  ffo/oi  (x,  y) 


hz 

Gz 


/(a:,2/)  +  j  j  ^  (a:,2/.a:o,2/o)  - (a:^2/^a:o,yo) 

D 


7-i(a:o,yo) 


+ 


^  (x,y,io,yo)  +  ^  {x',y',XQ,yQ) 


{xo,yo)  f  dxodyo  =  crofo2{x,y) 


(7) 


where 


foi{x,y)  =  Real  |  ^  -  \l -  <A 

-  pI'’''!'’  (4'*  -  V 4'**  - } 

,W4-'  (4"  -  V4^) 


/02(2:,  y)  =  Real  < 

I  Ml  -  P2 


-  92^ 


(4->-^ 


a^o  =  yo 

y'o  =  b-  XQ 

x'  =  y 
y'  =  b-x 


(8) 


The  parenthetical  numerical  superscripts  indicate  layer  number. 

As  in  reference  7,  the  kernels  Qy^  and  Rjk  {j,k  =  1,2),  which  have  logarithmic  singularities,  are  known 
and  are  integrable  in  the  domain  D.  Thus,  equations  (7)  can  be  treated  as  Fredholm’s  integral  equations  of  the 
second  kind  and  can  be  solved  numerically  by  reducing  them  to  a  system  of  linear  algebraic  equations.  If  the 
integral  signs  are  replaced  by  summations  and  collocation  is  used,  the  following  system  of  linear  equations  in 
the  unknowns  Tx{x,y)  and  Ty{x,y)  is  obtained: 

h  ^ 

-p^Txixj,  yj)  +  V  [5ii  {xj,  yj,  xo„,  yo„)  Tx  {xo„,yQj 

n=l 

+  Si2  (3:y, yy, xo„, yo„)  Ty  (xo„, yo„)]  AxnAy„  =  (tq/oi  {xjiyj) 

h  ^ 

-^Ty  (xy.yy)  +  V  [521  {xj,yj,X0„,y0n)  Tx  (xOn^yOn) 

n=l 

+  S22  {xj,  yj,  X0n,yQn)  Ty  Nn.  yOn)]  ^Xn^yn  =  CTo/02  {xj,  yj) 


{j  =  hN) 


(9) 


where 


Qii  ,  R22 


Sl2  = 


Qi2  _  R2I 

hi  h2 


(10a) 

(10b) 
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(10c) 


and  N  is  the  number  of  collocation  points, 
and  B,  respectively. 


o  Q21 

^^21  =  -7 - , 

hi  h2 


Ri2 

h2 

R\i 


Q  Q22  , 

*22  =  -r~-  +  , 
hi  h2 


(lOd) 


The  terms  Qji^  and  Rjj^  {j,  k  =  1,2)  are  defined  in  appendixes  A 


The  Stress-Intensity  Factor 

In  the  present  problem,  because  of  symmetries  in  geometry  and  loading,  the  shear  component  of  the  stress- 
intensity  factor  is  zero.  The  normal  component  is  found  by  combining  the  effects  of  the  crack-face  pressure 
and  the  distributed  body  forces  Xi{x,y)  and  Yi{x,y)  acting  on  layer  1.  The  stress-intensity  factor  may  be 
expressed  in  terms  of  the  unknowns  Tx  and  Ty  as  follows: 

fci  (ao)  =  o-Qi/a-l- y  [Ti{x,y,xo,yQ)Tx{xo,yo)  +  T2{x,y,xo,yo)Ty{xo,yQ)]dxodyo  (11) 

D 

where 


Ti  {x,y,X(i,yQ)  =  Fj  +  Fj 

T2  (2,2/,  XQ,  yo)  =  <7^2^  ^Fj  \  aj  +  Fj 


Fi{w,a)  = 


—  c?  —  [w  —  ao) 


w  ~  ao 


“)] + 42  [^/  (4'^ 


oj  +Fi 

«)  +  Fl 


and  00  is  replaced  by  -|-a  for  the  right  crack  tip  and  by  -a  for  the  left  crack  tip.  Complete  details  of  the 
derivation  of  the  equation  for  stress- intensity  factors  are  given  in  appendix  C. 

The  solution  of  equations  (9)  gives  the  shear-stress  distribution  in  the  adhesive  layer  Tj^x,y)  and  Ty{x,y). 
Using  these  adhesive  shear  stresses  in  equations  (11),  we  can  determine  the  stress-intensity  factors  at  either 
crack  tip.  The  numerical  solution  to  the  problem  is  discussed  in  the  next  section. 


Numerical  Solution 

A  key  item  in  the  analysis  is  the  method  of  in¬ 
tegration  used  for  the  system  of  integral  equations 
represented  by  equations  (7).  In  theory,  these  equa¬ 
tions  require  integration  over  the  infinite  domain. 
However,  because  of  the  complicated  nature  of  these 
equations,  a  closed-form  integration  is  probably  im¬ 
possible.  Consequently,  the  system  of  Fredholm’s  in¬ 
tegral  equations  is  solved  with  standard  numerical 
techniques.  This  is  done  by  dividing  the  domain  into 
cells,  assuming  the  unknown  functions  Tx  and  Ty  to 
be  uniform  in  each  cell,  and  then  using  a  numerical 
scheme  to  evaluate  integrals.  In  this  way  the  sys¬ 
tem  of  integral  equations  is  reduced  to  the  system  of 
algebraic  equations  given  by  equations  (9).  Two  im¬ 
portant  aspects  of  the  numerical  integration  are  the 
domain  of  integration  and  the  logarithmic  singulari¬ 
ties  in  the  kernels. 

In  order  to  carry  out  the  numerical  analysis,  the 
size  of  the  infinite  domain  must  be  restricted.  This 
restriction  is  possible  because  the  disturbance  caused 
by  the  presence  of  the  crack  decays  with  distance 


from  the  crack.  A  convergence  study  has  been  con¬ 
ducted  to  ensure  that  restricting  the  domain  does  not 
appreciably  affect  the  accuracy  of  the  results.  Thus, 
the  size  of  the  domain  D  is  determined  iteratively, 
starting  with  a  small,  coarse-mesh  domain  and  in¬ 
creasing  the  domain  extent  and  mesh  refinement  un¬ 
til  no  significant  changes  occur  in  the  quantities  of 
interest. 

The  second  aspect  of  the  integral  equations  re¬ 
quiring  special  attention  concerns  the  kernels  ^jk 
which  appear  in  equations  (9).  These  kernels  con¬ 
tain  logarithmic  singularities  which  must  be  evalu¬ 
ated  separately.  This  is  done  by  performing  an  exact 
integration  of  the  singular  portion  of  the  kernels  over 
the  segment  Aa:„.  In  effect,  this  integration  applies 
the  load  over  a  finite  length  rather  than  at  a  single 
point,  thereby  removing  the  singularity. 

Discussion  of  Results 

To  gain  insight  into  the  behavior  of  the  stringer 
panel,  the  analysis  was  performed  for  numerous 
combinations  of  material  and  geometric  parameters. 
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Three  different  materials  were  assumed  for  the  sheet. 
Their  properties  are  given  in  table  L  Material  4,  also 
shown  in  table  I,  was  used  to  model  the  stringer. 
These  materials  span  a  wide  range  of  properties  and 
are  representative  of  frequently  used  graphite/ epoxy 
laminates  (ref.  12). 

Figure  5  shows  the  normalized  stress- intensity 
factor  as  a  function  of  6/a,  the  ratio  of  the  stringer 
edge  distance  to  the  crack  half-length.  Here,  the 
crack  half-length  a  is  held  constant  and  6  is  increased. 
In  this  and  all  subsequent  figures,  the  stress-intensity 
factors  have  been  normalized  with  ao\/u,  the  stress- 
intensity  factor  for  a  cracked,  infinite  sheet  with  a 
uniform  pressure  cro  applied  to  the  crack  face.  As 
shown  in  this  figure,  once  the  right  crack  tip  is  no 
longer  beneath  the  stringer,  the  reduction  in  k\  is 
minimal.  Also,  the  effect  of  the  stringer  on  the 
left  crack  tip  is  much  less  than  that  on  the  right 
crack  tip.  This  is  to  be  expected,  since  the  left 
crack  tip  is  farther  from  the  edge  of  the  stringer. 
For  large  values  of  6/a,  both  the  right  and  the  left 
crack-tip  stress  intensities  approach  1.00  (the  value 
when  no  stringer  is  present),  again  as  expected. 
Changing  the  material  properties  of  the  sheet  does 
not  greatly  affect  the  value  of  fci,  although  the  lower 
the  longitudinal  stiffness  of  the  sheet  material,  the 
greater  the  reduction  in  the  stress  intensity  due  to 
the  presence  of  the  stringer. 

Figure  6  presents  data  similar  to  those  shown  in 
figure  5,  but  in  a  slightly  different  manner.  In  this 
case,  6  is  held  constant  while  a  is  increased.  The 
stress-intensity  factor  is  not  affected  by  the  presence 
of  the  stringer  until  the  crack  tip  is  very  close  to 
the  edge  of  the  stringer,  and  no  significant  reduction 
in  ki  occurs  until  the  right  crack  tip  is  beneath  the 
stringer  itself.  As  in  the  previous  figure,  the  stress- 
intensity  factor  at  the  left  crack  tip  is  not  significantly 
affected  by  the  stringer,  and  changing  the  material  of 
the  sheet  does  not  have  a  large  effect  on  fci. 

In  figures  7  and  8,  the  effects  of  varying  the  stiff¬ 
nesses  of  the  stringer  and  the  adhesive  are  examined. 
In  these  figures,  material  1  is  used  for  the  sheet  and 
material  4  for  the  stringer;  two  ratios  of  6/a,  1.2  and 
0.5,  are  used.  In  these  figures,  only  results  for  the 
right  crack  tip  are  shown,  since  any  effects  on  the 

left  crack-tip  stress-intensity  factor  are  minimal. 

(2) 

The  effect  of  varying  the  stringer  stiffness  Ey  ^/i2 
is  illustrated  in  figure  7.  The  normalized  stress- 
intensity  factor  for  the  right  crack  tip  is  plotted  as 
a  function  of  the  ratio  of  the  stringer  stiffness  to 
the  sheet  stiffness.  This  figure  shows  that  as  the 
stringer  stiffness  approaches  zero,  the  normalized  fci 
approaches  1.00,  the  solution  for  a  cracked  sheet  with 


no  stringer.  As  the  stringer  stiffness  increases,  the 
curves  approach  a  lower  bound  on  fci.  This  indicates 
that  past  a  certain  point,  no  further  reductions  in  ki 
can  be  obtained  by  increasing  the  stringer  stiffness. 
The  value  of  this  lower  bound  on  fci  is  a  function  of 
the  adhesive  stiffness  and  6/a.  Changing  the  stringer 
stiffness  has  much  less  effect  on  the  stress-intensity 
factor  for  the  case  of  6/a  =  1.2  than  for  6/a  =  0.5. 

The  effect  of  varying  the  adhesive  shear  stiffness 
Ga/hs  is  illustrated  in  figure  8.  The  normalized 
stress-intensity  factor  is  plotted  as  a  function  of  the 
ratio  of  the  adhesive  shear  stiffness  to  the  shear 
stiffness  of  the  sheet.  As  the  adhesive  shear  stiffness 
approaches  zero,  the  value  of  ki  tends  toward  1.00,  as 
expected.  As  the  adhesive  shear  stiffness  increases, 
k\  decreases.  For  6/a  ”  1.2  the  decrease  is  small, 
since  the  effect  of  the  stringer  is  negligible.  However, 
for  6/a  =  0.5,  the  decrease  in  fci  is  substantial. 
Although  these  points  lie  off  the  scale  shown  in 
figure  8,  for  very  large  values  of  G3  (on  the  order 
of  10^^  MPa),  the  stress- intensity  factor  does  not 
approach  zero  (as  might  be  expected)  but  rather 
some  nonzero  minimum  value.  This  lower  bound  is 
a  function  of  the  stringer  stiffness  and  6/a.  Figure  8 

also  shows  that  for  values  of  G^hi/G^yh^  greater 
than  about  4.0,  increasing  the  adhesive  shear  stiffness 
has  little  effect  on  the  stress- intensity  factor. 

Concluding  Remarks 

This  report  has  presented  the  formulation  of  an 
analysis  for  a  cracked  orthotropic  sheet  reinforced 
with  an  adhesively  bonded  orthotropic  semi-infinite 
sheet.  This  configuration  was  assumed  to  represent  a 
bonded  stringer  panel  when  the  crack  is  located  close 
to  the  edge  of  the  stringer.  The  adhesive  bond  was 
also  assumed  to  be  intact.  A  parametric  study  was 
conducted  to  determine  the  sensitivity  of  the  system 
to  material  properties  and  specimen  configuration. 
The  parameter  having  the  greatest  influence  on  the 
stress-intensity  factor  was  the  distance  from  the  near 
crack  tip  to  the  edge  of  the  stringer.  Unless  the  crack 
was  very  close  to  or  under  the  stringer,  the  reduction 
in  fci  was  minimal,  and  the  farther  the  crack  was 
under  the  stringer,  the  greater  the  reduction  in  fci. 
Increasing  the  stiffness  of  the  stringer  or  of  the  adhe¬ 
sive  also  resulted  in  a  reduction  in  the  stress- intensity 
factor. 


NASA  Langley  Research  Center 
Hampton,  VA  23665 
April  11,  1985 
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Appendix  A 

Derivation  of  Displacement  Functions  for  Infinite  Cracked  Orthotropic  Sheet 

This  appendix  presents  the  derivation  of  the  displacement  functions  for  the  infinite  cracked  orthotropic 
sheet  (layer  1).  If  we  follow  Lekhnitskii  (ref.  10)  and  ignore  an  arbitrary  rigid-body  displacement  vector,  the 
displacements  due  to  the  crack-face  pressure  cq  (shown  in  part  B  of  figure  3)  can  be  written  as  follows: 


(6)  Tj  ,  j  ^0 

‘‘  "  UW-„w 

_  Real  / 


(.!■>  - 


P2  >1 


zk  '  —  \  zX  —a 


(1)2 


(4^^  - 


(1)' 


(Al) 


where 


a  =  crack  half-length 

(y  =  i,2) 


Tp{k)  ) 


{k) 

a  ■  = 

^3 


.(fc) 


{k)  {k) 

"w  ■  ^ 


[j  =  1,  2;  {k)  =  Layer  number) 
{j  =  1,  2;  (/c)  =  Layer  number) 


The  terms  3,re  the  roots  of  the  characteristic  equation 


where  Ex  and  Ey  are  the  Young’s  moduli,  Gxy  is  the  shear  modulus,  and  i^xy  and  Uyx  are  the  Poisson’s  ratios 
of  the  orthotropic  sheet. 

To  determine  the  displacements  for  the  distributed  loading  shown  in  part  C  of  figure  3,  start  with  the 
displacements  for  a  concentrated  force  located  at  the  point  zq  =  xq  +it/o*  The  displacements  for  a  concentrated 
force  can  be  written  as 


where 


u%x,y)  =  2  Real  (4^^] 

v‘'{x,y)  =  2  Real  +4^42  (4’^)] 


(1) 

— (ly  +-BiCd 

“  ^2 

+  CiFd  {wi^\z\^\a)+DiFD 
-1-  Al  log  4-  Cl  log  -  Wg^^) 

(1) 

^2  (4^4  (1)  +BiFd 


Mi  -  M2 


(1)  Jl) 

4  ’^2 


CiFd  -|-T>iCd 

Bi  log  ~  '*"4^^) 


■»)] 


(A2) 
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F£)(ty,z,a)  -  log  (^\/w^ 


a^y  —  a? 


-\~wz  —  j  —  log  ^\/ 


(1)  .  (1) 
w\^  =  XQ-h  'yo 

(1)  (1) 

wy  =  XQ-  III  'yo 

Ai  =  C^^Xi+C^^Yi 

Cl = ci;)xi  -  cJjVi 


(1)  ,  (1) 

w\’  =XQ  +  ’yo 

(1)  (1) 
wi  =xo-  tJ-2  yo 

Bi  =  cWXi+cgVi 

Di  =  c(;)Xi  -  cgVi 


.(i)_  1 

(^1  +  4^42^ ^2^4  4^^  +  4'^^] 

11  “  ^  •  / 
ZTTZ  {fl 

.(i)_  \ 

W-4«) 

(4'^^ + 4'^4  4^^ + 4'^^4'^^ + 4y  ] 

12  ~  „  .  ( 

2'ki  I /:x 

.(i)_  ^2  L 

4->)(i+4V1'¥.'')  +  p?>+4'>] 

21  “  „  .  / 
2711  \fX 

.0)_  ^2  1 

[((i«+,.«)4'>+4¥f +45] 

Details  on  the  derivation  of  <f>l  and  (j)^  (4^^)  found  in  reference  13. 

Since  the  body  forces  Xi  and  Fi  are  continuous  functions  of  zq  in  the  domain  D,  the  displacements  and 
v^'  may  be  used  as  Green’s  functions  to  determine  ^  and  \  The  displacements  ^  and  Uj  ^  due  to  the 
body  forces  Xi  and  Yi  distributed  over  the  domain  can  be  written  as  follows: 


\\x,y)  —  I  {x,y,xo,yo)  dxQdyo 


■  j  j  [Qii  (a:,  y,  XQ,  yo)  Xi  {xq,  yo)  +  Qn  y,  ^o,  yo)  Fi  (xo,  yo)]  dxo  dyo 


’{x,y)—  /  /  v‘^(x,y,  xo,yo)  dxodyo 


=  j  j  IQ^I  {x,y,xo,yo)Xi{xo,yo)  +  Q22  (x,y,a:o,yo)Fi  (xo,yo)]  dxodyo 


where  Qn,Qi2,  Q2I)  Q22  are  found  by  substituting  equations  (A2)  into  equations  (A5)  and  then  combining 
terms  and  simplifying.  These  expressions  are  shown  below. 


9 


gn  (x,.,xo,.o)  =  2  Real  (.f .  a)  +  C^F, 

.  f^l  i^2 


-  («>i”.4'>,-) +41Vc  H‘'.4‘>.«) 

Mi  -  M2 

+  c<'>oH‘>,4''.«)+4i’foK’.4‘>.«)l 

2M9  P.  R..  +  ^21  j  + 


-  4^^ 


24»pWfc!;)+c« 


4**  -  4^' 


log  (  Y  4^^  -  a2  + 


+  Pl^^^ll^  [log  +  log  {z[ 

+  [log  (4^^  -  ^2^^)  +  log  (4 


(1)_  Jl) 

1  ^3 


.(i)„(i) 


gi2  (..  y,  .0,  .0)  =  2  Real  [cglFG  (.14  a)  +  c4^ 

I  Ml  M2 


^Vi^G  (w4>4^^®)  -^22^G 


.,(1)  ,(1) 


.(1)„(1) 


-^^[c4)FG(.i444«)  +  c44G(. 

lSFal4\^i^\a)-cSFo{J^\z?,a)] 


(1)  (1) 
W2  ,^2 


+  ^4^2^  [log  (24  -  w;4)  -  log  (^4  -  ■^4)] 

[log  (4^^  -  -  log  (4^^  -  I 


,(l)..(i) 


Q21  {x,y,xo,yo)  =  2  Reali  ^^i)  [^11^ (w^, a)  +  cI^IFg  («;4,a4,a) 

I  Ml  “"M2 


+  cJilFG  4'^  «)  +  ( 


w^^\z[^\a 


Ml  -  M2 

+  C^Fg  {wi^\zi^\a)] 

o„(l)Jl) /^(l)  ,  f7(l)\  r  // _ 

2m2  9i  [Cn  +C21) 


21  ^G  (w4.4^^«) 


m4  “  m4 


.4^  -  a2  +  41) 


(Equation  continued  on  next  page.) 


Appendix  B 


Derivation  of  Displacement  Functions  for  Semi-Infinite  Orthotropic  Sheet 

This  appendix  presents  the  derivation  of  the  displacement  functions  for  the  semi-infinite  orthotropic  sheet 
(layer  2).  The  only  forces  acting  on  the  semi-infinite  sheet  are  the  distributed  body  forces  X2  and  Y2,  as  shown 
in  part  C  of  figure  2.  In  a  procedure  similar  to  that  used  in  reference  11,  displacements  due  to  a  pair  of  body 
forces  X2  and  located  at  point  ^0  ~  ^0  oriented  as  shown  in  figure  4,  can  be  written  as  follows: 


where 


where 


{x',y\xQ,yQ)  =  2 Real  (4^^ 

(i',y',xo,yo)  =  2 Real  (4^^ 


(2)  /  J2) 


A2  log  ^2  log 


^2^^  - 


{  -  4^4  [^2  log  -  wf 


+  G2  log  -  4^41  +  (^2^^  “  4^^ 


X  [S2  log  [z^-^  -  + D2  log  I 

<^2^^  (4^4  ^  (^2^^  ~  ”*^2^0  '''  ~  “’4^0 

+  jiT  m  { -  "I"’)  h  i»s  (4"’  - 

Ml  M2 

+ C2  log  {zp  -  4]  +  (^2^  -  mi^4 

X  [r2  log  (z^^^  -  -L  n.  f  ,(2)  _  .7,(2)M  \ 


^2  + 


D2  log  41 


=  ^0  +  ^^f'‘yo 


(i  =  i,2) 


=  Xq  +  4o 

4^^  =  -^0 + 4o 
A2  =  c[fx'2  +  cgV2' 
C2  =  -c[fx'o  +  dfY^ 


,,,(2)  -  J  >  ,,(2)  / 

^2  “  ^0  +  P2  2/0 

4^^  =  -2:0  +  P2^4o 

r2  =  4?x^+4?^2 
7?2  =  -4l^X4  +  cgV2' 


and  42^’  44  42^  defined  in  equations  (A4).  Details  on  the  derivation  of 

4>^2^  (4^4  found  in  reference  13. 

As  before,  if  the  body  forces  X2  and  Y2  are  continuous  functions  of  Zq  in  the  domain,  the  displacements 
and  1;^,  given  by  equations  (Bl),  may  be  used  as  Green’s  functions  for  the  displacements  and  due  to  the 
distributed  forces.  Thus,  the  displacements  can  be  found  as  follows: 


4  {x',y')  =  yy  w*  {x',y',XQ,yo)  dx’ody'o 
D 

~  I  I  "*” 

D 

4  (x',y')  =yy (a:',  2/',  4 >2/0) 

D 

=  j j  [i?2i  (a:',  y',  4,  y'o)  ^2  (xq,  yo)  +  R22  {x',  y',  Xq,  y'o)  Y2  (xg,  yg)]  dxg  dyg 


(B3) 


where  iin,  i2i2,  i22l7  i222  are  found  by  substituting  equations  (Bl)  and  (B2)  into  equations  (B3)  and 

then  combining  terms.  The  complete  expressions  for  i?ii,  i2i2,  i22l5  and  R22  are  given  below. 

Equations  (B3)  give  the  displacements  of  the  semi-infinite  sheet  for  the  loading  and  the  coordinate  system 
shown  in  figure  Bl(a).  However,  the  displacements  U2  and  V2  required  for  the  solution  are  defined  with  respect 
to  the  coordinate  system  and  are  due  to  the  forces  X2  and  I2  shown  in  figure  Bl(b).  Thus,  from  figure  Bl, 


where 


U2{x,y)  =  -4  (x',y') 


U2(x,y)  =  4  {x',y') 


x'  =  y  y'  —  h  —  X 

X'2  =  Y2  Y^  =  -X2 


(B4) 


The  displacements  in  the  semi-infinite  sheet  (layer  2)  are  now  found  by  combining  equations  (Bl)  and  (B3) 
and  substituting  into  equations  (B4).  Note  that  the  material  properties  of  layer  2  are  defined  with  respect  to 
the  x',y'  coordinate  system  (fig.  4). 

Shown  below  are  the  complete  expressions  for  the  functions  i^n,  i2i2,  i?2l!  ^22)  which  appear  in 

equations  (B3).  Further  details  on  the  formulation  of  these  expressions  can  be  found  in  reference  13. 


i?ii  [x',  y',  x'g,  y'o)  ^  2  Real  (pf  [log  ^  -  log  {zf^  -  wf 

+  [log  (4^^  “  ^2^)  ~  log  (4^^  “  ^)] 

+  T#U5  { -  4^’)  <5!?  h  -  ■'’*  -  '*“)] 

f^2 

+  (4^1  -  ,?>)  cp)  [log  (.f)  -  )  -  log  (4«  -  )] } 

+  T5f^  {(-?’  -  d^’)  <?!?  h  (4’’  -  d^’)l 

Pi  -  M2 

+  -  Pi^^)  1^2?  [log  (4^^  -  '^2^)  -  log  (4^^  -  41 }) 
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Ri2  {x,  y',  x'q,  y'o)  =  2  Real  log  ^  +  log  [zf'^  -  wf>"^ 
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P2  ^1 

(m2^^  -  4^4  ^2?  [log  ^  -  4^4  (^1^^  “  ^4^41  } 

-#py  {(pf>  -  P) <?!!>  h  P'  - ’=’®)  +  - ’^"’)] 

Ml  Mo 


^^2 

4^^  -  /xf  4  1^22^  [log  (4^^  -  "^2^4  (4^^  “  ^4^^ 
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Appendix  C 

Derivation  of  Stress-Intensity  Factor  Equations 

This  appendix  presents  the  derivation  of  the  stress-intensity  factor  equations.  In  the  present  problem, 
because  of  symmetries  in  geometry  and  loading,  the  shear  component  of  the  stress- intensity  factor  is  zero. 
The  normal  component  is  found  by  combining  the  effects  of  the  crack-face  pressure  (tq  and  the  distributed 
body  forces  Xi{x,y)  and  Y\{x,y)  acting  on  layer  1.  For  an  orthotropic  material,  the  normal  component  of  the 
stress-intensity  factor  may  be  written  as  follows  (ref.  14): 

fci  =  2n/2^^  lim  \] -  a  4>'i  (4^0  dxodyo  (Cl) 

n 


fci  =  2n/2^^  lim  +  // <^1"  (4^0  dxodyo 

y-2  L  V 

The  function  >  the  stress  potential  due  to  the  crack- face  pressure  (tq,  is  given  by  the  following: 

,,  /  (i)\  ^  / 1  _  4^^ 


4^'  -  o2 


The  function  stress  potential  (derived  in  ref.  13)  due  to  the  body  forces  Xi  and  Yy,  is  given  by 

the  followinff: 


’^‘'("F)  =  (■)""  (1)  (T’.-i".") 

H  ~  ^2  \j  z\  ^ 

+  BiFh  +CiFh 

+  DiF„  (z[  \w[  \a)j  +  ^(1)  (1)  +  (1)  _  (1) 


where 


Fjj{z,w,a) 


_  y  —  0?  —  [z  —  w) 


^  n  v'^5  /V 

[z  —  W) 

and  Ai,  Bi,  Ci,  and  Di  are  defined  in  equations  (A3).  Substituting  equations  (C2)  and  (C3)  into  (Cl),  taking 
the  limit,  and  combining  terms  yield  the  following  stress-intensity  factor: 

ki=(7o\/d-  ^  j  j  [Ti  (x,  y,  XQ,  yo)  Xi  (xq,  yo)  +  T2  (x,  y,  xq,  yo)  Yi  (xq,  yo)]  dxQ  dyo  (C4) 


where 


Ti  (x,  y,  XQ,  yo)  -  ^Fj  -I-  Fj  aj  -I-  ^Fj  -f  Fj  (^4^^  a) 

T2  [x,  y,  XQ,  yo)  =  C^2  [-^1  «)  +  ^Fj  aj  -I-  Fj  (^^4^^  oj 


Fj{w,a)  = 


I  vfl  —  a2  —  w  +  a 


Through  the  use  of  equations  (2),  equation  (C4)  may  be  written  in  terms  of  the  unknowns  r^and  Ty  as  follows: 

fci  =  cro-/a  +  j  j  [Ti  (x,  y,  xq,  yo)  Tx  (xo,  yo)  +  Tz  (x,  y,  xo,  yo)  Ty  (xo,  yo)]  da:o  dyo  (C5) 
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Since  the  loading  on  layer  1  (as  shown  in  fig.  3)  is  not  symmetric  about  the  y-axis,  equation  (C5)  is  valid  only 
as  zi  approaches  a,  that  is,  at  the  right  crack  tip.  To  find  the  stress- intensity  factor  at  the  left  crack  tip,  the 
function  (j)^  found  for  the  forces  Xi  and  Yi  located  in  the  left  half-plane  (>2^0  —  —xq  -\-iyo)  • 

With  the  new  form  of  (j)^  same  procedure  as  above  is  used  to  find  ki  at  the  left  crack  tip.  The 

function  (j>[  does  not  change,  since  it  does  not  depend  on  the  location  of  the  concentrated  forces  Xi  and 

Yi.  The  stress-intensity  factor  at  the  left  crack  tip  is  thus  found  to  be 

y,  ^0,  yo)  ‘Tx  (3:0,  yo)  +  {x,  y,  xq,  yo)  Ty  (xq,  yo)]  dxQ  dyo  (C6) 


where 


T3  {x,  y,  xo,  yo)  =  +  Fjr  (w^^\  aj  j  +  ^Fp  aj  +  Fp 

T4  (x,  y,  xo,  yo)  -  [i^F  «)  +  «)]  +  cS  [^>  (4'^  «)  +  ^f  «)’ 


Fp{w,a) 


IvP'  —  0?  —  (ly  +  a) 
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TABLE  I.  MATERIAL  PROPERTIES 


Material 

Ex, 

GPa 

Ey, 

GPa 

^xy 

Gxy, 

GPa 

Graphite  /  epoxy 
laminate 

1 

51.40 

51.40 

0.3065 

19.67 

[0/  ±  45/90]s 

2 

19.75 

19.75 

.7336 

33.65 

[±45]2s 

3 

129.40 

10.86 

.3118 

5.70 

[90] 

4 

10.86 

129.40 

.0262 

5.70 

[0] 

A. Total  loading 


B.  Crack  pressure 


C.  Distributed  loading 


Figure  3.  Superposition  for  sheet  (layer  1). 


Figure  5.  Normalized  stress-intensity  factor  as  function  of  edge  distance,  a  —  10  mm;  G3  =  460  MPa; 
h2/hi  =  0.5;  h^/hi  =  0.005. 


0  12  3 


a/b 

Figure  6.  Normalized  stress-intensity  factor  as  function  of  crack  length.  6  =  10  mnrGs  =  460  MPa; 


b/a=  1.2 


kj{+a) 

0^ 


b/a  =  0.5 


3.0  4.0  5.0  6.0  7.0 

P^2)h 

y  1 


Figure  7.  Normalized  stress- intensity  factor  as  function  of  stringer  modulus. 


b/a  =  1.2 


kj{+a) 


4.0  6.0 


b/a  =  0.5 


10.0  12.0 


h,  G  ' 

3  xy 


Figure  8.  Normalized  stress- intensity  factor  as  function  of  adhesive  shear  modulus. 
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